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ABSTRACT 

Context. The detection of an extra component in GRB 060904B X-ray spectra in addition to the standard single power-law behaviour 
has recently been reported in the literature. This component can be fit with different models; in particular the addition of a spectral 
line provides the best representation. 

Aims. In this paper we investigate the physical properties that the surrounding medium must have in order to produce a spectral feature 
that can explain the detected emission. 

Methods. We analyse and discuss how and if the detected spectral excess fits in different theoretical models developed to explain 
the nature of line emission during the afterglow phase of Gamma-Ray Bursts (GRBs). Trasmission and reflection models have been 
considered. 

Results. Given the high value (» 1) of the Thomson optical depth, the emission is likely to arise in a reflection scenario. Within 
reflection models, the external reflection geometry fails to predict the observed luminosity. On the contrary, the detected feature can 
be explained in a funnel scenario with typical opening angle 9 ~ 5°, Nickel mass ~ 0.1 M Q and T = 10 6 K. For 9 ~ 20°, assuming the 
reprocessing material to be the SN shell, the detected emission implies a Nickel mass ~ 0.4 M G at T ~ 10 7 K and a metallicity ~ 10 
times the solar value. If the giant X-ray flare that dominates the early XRT light curve is identified with the ionizing source, the SN 
expansion began ~ 3000 s before the GRB event . 

Key words. GRB: X-Ray afterglow 



1. Introduction 

A direct observation of Gamma-Ray Burst (GRB) central en- 
gines is not possible; however, it is possible to infer something 
about the nature of their progenitors if indirect probes of the 
physical conditions, structure and composition of the material 
at the GRB site can be found. X-ray spectroscopy is one of such 
probes. 

Up to the present time, the overwhelming majority of X-ray 
spectra of afterglo ws detected by Swif t shows a power law be- 
haviour, (see e.g. O'B rien et al.l l2006) result of a non thermal 
emis sion, the le ading candidate for which is syncrotron emis- 
sion (Piran 2005 and references therein). Any deviation from this 
standard behaviour would be extremely interesting. 

Detections of emission and absorption features in addition 
to the basic absorbe d power law have been claimed in a number 
of observations (see ISako et al.ll2005l their Table 2 for a com- 
prehensive summary). Most are interpreted as Fe Ka emission. 
The emission line radiation is supposed to be produced through 
reprocessing of the burst radiation: for this reason its luminosity, 
spectral shape and wavelength shift carry important information 
not only about the geometry and structure of the surrounding 
medium, but also about the ph ysical conditions at th e GRB site 
and the same GRB energetics (Ghisellini et al. 2002). However, 
the statistical significance of th e claimed detections has been 
questioned by Sak o et all (120051) : a large fraction of the claims 
are based on using the F-test, a statistical method that gives both 
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false positives and f alse negatives when ap plied to emission or 
absorption features dProtassov et alj|2002l) . After a re-analysis 
of the data based on Monte Carlo simulations lSako et alJ (2005) 
ruled out all of the reported features. 

Recently, new d etec tions have been claimed bv lButleri (120071) 
and iMoretti et aT] d2007l) (hereafter M07). In particular, M07 
focussed on spectra of GRBs with know redshift detected by 
Swift in the period April 2005-January 2007 with very steep 
spectra (photon index F > 3, value obtained with an absorbed 
single power law -hereafter SPL- fit). Out of 13 spectra with 
more than 2000 photons satisfying this condition, highly signifi- 
cant deviations from the SPL spectral model have been found in 
4 cases. These spectra belong to GRB 060502A, GRB 060729, 
GRB 060904B and GRB 061 1 10A. These bursts have been de- 
tected by Swift. In one case (GRB 060904B) the excess can be 
modeled by the emission line of highly ionized Nickel (see Fig. 
ID- 

We refer the reader to M07 for details about the sample se- 
lection and the statistical methods used. 

In this paper we analyse and discuss how and if the spec- 
tral excess detected in GRB 060904B fits in the different models 
developed to explain the production of line emission during the 
afterglow phase of GRBs. In particular, the theoretical implica- 
tions of the possible detection of Nickel emission » 200 s (see 
Fig. |2]) after the onset of the GRB event will be explored. 

This paper is organised as follows: the spectroscopic iden- 
tification of the line emission is discussed in Sect. [2] In the 
same Sect, we derive simple constraints to the theoretical mod- 
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Fig. 1. Upper panel: GRB060904B spectrum observed in the rest frame time interval 138 - 185 s. Solid thick black line: best 
fit absorbed simple power-law model (SPL): photon index F = [jg; neutral absorber column density at the source redshift 

N H ,z = (0.69!° {J;}) x 10 22 cm" 2 . Thin coloured line: best fit absorbed power-law plus gaussian model (GAU): F = 3.67+° 0°; Mt,z = 
(0.74+° °^) x 10 22 crrT 2 ; the values of the parameters related to t he gaussian component are r eported in Tab.Q] The Galactic hydrogen 
column density A^h.mw has been fixed to the value reported in lDickey & Lockmanl (1 19901) along the GRB direction (M07). Lower 
panel: ratio between the observed values and the model predicted values. (In color in the online edition.) 
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Fig. 2. Luminosity curve of GRB 060904B in the observed 0.3 - 
lOkeV energy band (XRT observations). Rest frame times as 
measured from the BAT trigger have been used. Grey dashed 
lines: interval of time during which the anomalous emission 
arises. 



els starting from purely geometrical considerations. We discuss 
the relevant physical quantities linked to the detected emission 
in Sect. [3] Theoretical models are derived in Sect. [4] We address 
the problem of the ionizing source in Sect. [5] while in Sect. [6] 
we derive limits to the GRB total energy reservoir starting from 
the detected line luminosity. Finally, our results are discussed in 
Sect. [7] Conclusions are drawn in Sect. [8] 



The data reduction and temporal analysis have been per- 
formed using the standard HEADAS software, v6.1, while for 
the spectral analysis we used XSPEC (vll.3). Unless other- 
wise stated, we quote errors at the 68% of confidence level. 
Through out this paper standard concordance cosmology is used: 
Q. m = 0.27, S1 A = 0.73, h = 0.7 1 . 



2. Spectroscopic identification and geometrical 
constraints 

For GRB 060904B, the statistically significant deviation from 
the SPL spectral model is best modeled by an additional 
Gaussian component (see Fig. [T}. After a Monte Carlo analysis, 
M07 report a multitrial significance of the excess greater that 
99.9%. We refer the reader to M07 for details. Table|T]shows the 
physical quantities related to the detected excess. Unless other- 
wise state d, rest frame values w ill be used (the burst is located at 
z = 0.703 : lFugazza et al.ll2006l) . 

The detected excess is centered at 7.85 keV in agreement 
with what we expect from highly ionised Nickel line emission: 
the k a emission of H-like (He-like) Nickel is expected to lie at 
8.10keV (7.81 keV). Such detection differs from previous claims 
of line emission in GRB afterglows because of the much larger 
lumi nosity involved (about 3 orders of magnitude higher; see 
e.g. lBottcherl f2004 and references therein), and because of the 
really short delay between the GRB event and the line detection 
(» 10 2 s). In previous claims of line detection the event has been 
reported to occur at later times, i.e. 10 4 - 10 5 s after the GRB 
detection. 



R. Margutti et al.: Anomalous X-Ray emission in GRB 060904B: a Nickel line? 



3 



Table 1. Best fit parameters and physical quantities related to 
the spectral excess detected in GRB 060904B X-ray afterglow. 
The spectrum has been modeled with a SPL with superimposed 
a Gaussian component to fit the line profile: time interval be- 
tween the GRB explosion (BAT trigger) and the detection of the 
spectral excess; observed duration of the anomalous emission; 
Gaussian width; Gaussian central energy; luminosity of the ad- 
ditional spectral component; equivalent width. Rest frame values 
are provided. (From M07). 



the emitting Nickel mass is given by M^i = 56m p NAt « 20 M . 
Q However, Nickel masses as large as « 20 M are difficult to 
account for in any reasonable physical scenario. We therefore 
conclude that each Nickel atom recombines more than once in 
At « 47 s, or f rec < 47 s. Using the expression of the recombina- 
tio n time of an hydrogenic ion of atomic number Z = 28 given 
bv lVerner & Ferlandl (1 19961) the previous expression turns into a 
condition on the electron density and temperature of the plasma. 
The allowed parameter space is shown in Fig. [3] 
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" Lower limit to the real duration of the emission because of a gap of 
observation between ~ 200 s and ~ 2000 s (rest frame values) after the 
GRB onset. Hence we can write 47 s < At < 2000 s. 



The link between GRB events and supernova (SN) explo- 
sio ns has been established in a number of observations (see 
e.g.lGalama et al.ll998tlStanek et al.l2003HMalesani et all2004 
Cam pana et al.ll2006l) . If we require the line emitting material to 
be the SN shell, then from the failed detection of Co features at a 
3cr confidence level we derive an upper limit on the remnant age: 
f remn ^ 6 days. ( The time evolution of Ni, C o and Fe abundances 
as modeled by Woosl ey"& Weaver! 1 19951 has been assumed). 
Adopting v s heii ~ 10 9 cms _1 as ty pical expansion velocity of 
the SN shell at early stages (see e.g. lPatat et al.l l2001). the maxi- 
mum distance travelled by the reprocessing material is therefore 
^max ~ 10 15 cm. Assuming a contemporaneous SN-GRB explo- 
sion this distance shrinks to « 4 x 10 12 cm. On the other hand, 
from purely geometrical considerations the detection of the line 
emission at fd et requires the material to be located within a dis- 
tance/? = cf det (l -cos0) _I a 4xl0 12 (f det /138 s)(l -coscpy 1 cm. 
Note that R > 10 16 cm only for <p < 1°, where <p is the angle be- 
tween the line emitting material and the line of sight at the GRB 
site. For these reasons we consider unlikely the possibility that 
the detected line emission originates in a distant (R > 10 16 cm) 
reprocessor scenario and we prefer nearby {R < 10° cm) repro- 
cessor models. In such models the duration of the line emission 
is primarily linked to the time for which the central engine is 
active. Unfortunately we have only a poor lower limit to this pa- 
rameter (At > 47 s) so that also geometrical factors can give 
important contributions to the observed line duration. (In the 
first ss 130 s the high level of the continuum emission prevents 
us from detecting any spectral feature with the same luminos- 
ity even if present, while no data have been collected by Swift 
between ~ 200 s and ~ 2000 s after the onset of GRB 060904B). 



3. Physical Conditions 

Emission Time Scale- A line photon is produced each time an 
electron recombines. The time scale of the emission process is 
therefore set by the slowest between the ionization and the re- 
combination processes. Even under the assumption of a ioniz- 
ing luminosity similar to the luminosity of the observed excess, 
it can be shown that f; on <k f rec , so that the ionization is es- 
sentially istantaneous. The observed line luminosity implies an 
emission rate of Nickel line photons of N — Liine/fNi ~ 10 55 s _1 
(ejsH ~ 8 keV); if each Nickel atom recombines only once, then 
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Fig. 3. Coloured area: electron number density-temperature pa- 
rameter space valid for f rec < 47 s for GRB 060904B. The dashed 
lines mark the 0.1 M Q < M Ni < 20 M Q region, where M Ni 
is the emitting Nickel mass. The shaded area represents the 
T < 3.2 x 10 5 K region. The electron temperature is likel y to be 
highe r than this value at this stage of the SN explosion (lArnettl 
1996). 



Thomson Optical Depth- The Thomson optical depth of the 
SN shell of radius R and total mass M is given by: 



1 cr T / M\ 



p e m v AnR 2 \ X 



X ) 



(1) 



where p e is the mean molecular weight; X = M N ;« /M; is 
the total Ni mass. The photoionization optical depth for He-like 
or H-like Nickel can be written as: 



g cr Ni Xfl e 

crj 56 



Tj 



(2) 



where £ is the ionised to total Nickel abundance ratio and cr Ni 
is the He-like or H-like Nickel photoioniza tion cross section 
(« 10~ 20 cirT 2 , from lVerner & Yakovlevl 19951) . It is easy to show 
that t t » 1 for R ~ 10 12 - 10 13 cm and M > 10 5 M , where 
M is the total mass of the remnant shell. While an efficient use 
of the H-like or the He-like Nickel ions requires tni > 1, a 
Thompson optical depth larger than one is a problem: if photons 
diffuse through a plasma cloud in which Thomson scattering is 
the predominant mechanism and if the cloud has a central photon 
source, then the average number o f scatterings experienced by a 
photon will be « 1/2 t 2 (see e.g. ISunvaev & Titarchukl [1980). 
With tt » 1 , the line would be smeared in the continuum and 
no detection would be possible. We must therefore require the 
line photons not to have crossed such regions along our line 
of sight. This in turn implies that the electron clouds are not 
isotropically distributed around the GRB location, in order to al- 
low the line photons to freely escape. A tt » 1 favours models 



1 Note that in this paper we refer to 56 Ni since this is the Ni isotope 
mostly produced during SN explosions. 
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in which the line is produced by reflection; in this case line pho- 
tons come from the thin layer of material with tn; ~ several or 
tt < 1 . Trasmission models instead require the emitting material 
to globally have tt < 1. 

Line Width- In this paragraph we consider the different phys- 
ical processes that could lead to the observed line broadening 
(w 0.5 keV, Table [TJ. A purely thermal broadening requires 
T ~ 10 12 K, a temperature that is difficult to account for in 
any reasonable physical scenario. At the same time this would 
also cause the recombination efficiency to drop. On the other 
hand, if the bulk expansion of the SN remnant is the source of 
the observed broadening, then a velocity of the emitting atoms 
w 25 000 km s _1 is required. It is notable that such velocities are 
indeed observed during the Arts stages of Hy pernova (HN) ex - 
plosions associated with GRB events (see e.g. lPatat et al.l l2001). 

An alternative broadening mechanism is the Compton scat- 
tering by free electrons of the emitting material. In high tempera- 
ture gas the first scattering will pro duce an average energy broad - 
ening AE/E « (2k b T /m e c 2 ) 1/2 dSunvaev & Titarchukl fl980l) . 
For the emission feature detected in GRB 060904B spectra, this 
directly translates into T * 10 7 K. It is however possible that the 
real electron temperature is lower and that line photons undergo 
numerous scatterings before reaching the observer. In this case 
the centroid of the line would be redshifted (line photons are 
mostly backscattered by colder electrons) with an average en- 

1 /2 

ergy shift given by AE » 0.12 A^ sc keV, where the N sc stands fo r 
the number of suffered scatterings dSunvaev & Titarchuk| [T980). 
At the same time the multiple scattering process would further 
smear the line out, making it more difficult to detect. It is there- 
fore reasonable that N sc ~a few. 

Given the poor statistics, it is not possible to discuss in detail 
the line profile. In the following, results will be discussed for 
electron temperatures in the range 10 6 - 10 8 K. 

4. The models 

In this section we discuss models in which the line is due to 
reflection (see Sect. [3). The underlying physical mechanism is 

fast ionizations and recom binations. 

Reflection models (see Bottcher 2004 for a review) require 
the presence of a dense optically thick medium surrounding the 
burst site. Moreover, an anisotropic distribution of the material 
is required in order to have a clean path where line photons can 
propagate and reach the observer without being smeared in the 
continuum (see Fig. 0). The ionizing flux is efficiently repro- 
cessed in line photons in a superficial layer of the material with 
tt ~ 1 (in order to avoid excessive Compton broadening) and 
TNi ^several, to efficiently reprocess the continuum radiation. 
The volume V em effectively contributing to the observed line 
emission is therefore given by V em ~ S AR or V em * S /o- T n e . 
The line luminosity can be expressed as: 

r «NiV em e N i " Ni5 e Niffi-(Z, T) 

Mine - - (J) 

free O" j 

where: 

o" T : Thomson cross section; 
n e : Electron number density; 

«Ni : Ni number density; 
S : Emitting surface; 

6Ni : Energy of a single line photon (« 8 keV); 

a r {Z, T) : Recombination coefficient as a function of the 
atomic number Z a nd temperature T (from 
IVerner&Ferlandlll996l) . [or r ] = cm 3 /s; 



In the following two alternative geometries of emission will be 
considered: an e xternal reflection geometry and a stratificated 
funnel geometry dVietri et al.ll2.00lh . 

4. 1 . External Reflection Model 

According to the external reflection model, the line emitting ma- 
terial is back-illuminated by burst or afterglow photons scattered 
by preburst stellar wind (we refer the reader to lVietri et alj 2001 
and references therein for details). In this model the line emitting 
material has been ejected at subrelativistic speeds along the pro- 
genitor equator in a simultaneous SN-GRB explosion. The lu- 
minosity L scatt back-scattered by electrons mov ing outward from 
the source is of the order of ( Vietri et al. 200l]): 

L scatt « 1.8 xlO 45 ^^ ergs" 1 (4) 

V w ,7 

where m Wj _5 = ffZw/(lO~ 5 M0yr _1 ) and the symbol m w stands for 
the progenitor mass loss rate; v w j = v w /(10 7 cms _1 ), being v w 
the wind velocity. It is notable that this is already 2 orders of 
magnitude lower than the observed line luminosity and « 4 or- 
ders of magnitude l ower than the expected ionizing luminosity: 
lLazzati et al.1 (12002) show that in optimal conditions a reprocess- 
ing efficiency of w 1% at most can be reached, so that for the 
emission feature detected in GRB 060904B L ion > 10 49 erg s" 1 
is required. A fine-tuning of the progenitor wind velocity and of 
the mass-loss rate is therefore required to explain the observed 
line luminosity. For this reason we consider this scenario un- 
likely. 




Fig. 4. Sketch of the geometry assumed for the reflecting mate- 
rial for the funnel model discussed in the text. 



4.2. Funnel Model 

In this model the line emitting material is provided by t he walls 
of a funnel excavated in a young SN shell dLazzatil2000l) . A cone 
geometry of opening angle 9 is assumed for simplicity (see Fig. 
@). Considering a power-law electron number density profile: 

«e = «,(£) " (5) 
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equation @ translates into: 



^line — 



56(T X 



fNi 2tt a T (Z, T) sin(0)niR 



Jr. 



r'^dr 



where: 



56 56 M 
The totaQNi mass can be written as: 

M* Ni = X M = X ^47ryL( e m p cos(6')n I R; ! ' J r 2- " dr) 



(8) 



where m p is the proton mass. This allows us to derive the line 
luminosity independently of the quantity n\R": 



Loot = M N ——— tan(0)- 1 

112cr T m p f °r 2 - ff dr 



(9) 



The previous relation assumes f rec > t lon . A rough estimate of 
the ionization time scale f; on for material located at a distance R 
from the photon source of luminosity L\ on is given by: 



iionC Ni 



(10) 



where e; on is the energy of the ionizing photons; the parameter Q 
allows for the possibility of collimated flux of ionizing photons. 
The recombination time scale can be written as: 



tree — 



1 



n e a r {Z, T) 



(11) 



For the funnel cone geometry, the condition f rec > f; on therefore 
becomes: 



^ion°~ Ni 



1 



fio„ ae lm aAZ,T) ni R"r 2 - a 



> 1 



(12) 



As before, inserting n\R" from equation ([8]l and assuming the 
collimation angle of the ionizing photons to be ~ 6 we finally 
derive: 



free (honO~m 2cOS(6») ft/ftp 1 



free _ / Mon°"f 



eionar (1 -cos(6>)) M r 2 



> 1 (13) 



The system of equations that defines the funnel model is there- 
fore: 



^e=^ Ni ^g)tan (( J)|^ 

^ion^Ni 2cos(fl) AV« P i rR„ 2 -a j . i 
6 iu „ff r (l-cos(fl)) M r 1 -" JR, 

R -Ri~ cAt 



(14) 



where in the last equation the duration At of the line emission 
is required to be comparable with the time needed by the source 
photons to cross the SN shell. The observed or fixed quantities 
are the following: 



2 The fraction of the total Ni mass actually contributing to the de- 
tected emission is indicated with M^i- 
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Fig. 5. Upper limits to the total mass of the shell derived from 
the second relation of the system [14] evaluated at r — R Q as a 
function of the funnel opening ang le. Solid lines: T = 10 6 K; 
dashed lines: T = 10 7 K; dotted lines: T = 10 8 K. Thick lines: 
a = 0; thin lines: a = 1. 



Lii n e : Line luminosity (1.10 x 10 47 erg s _1 ); 

e; on : Energy of the ionizing photons (the ionization 
threshold is ~ lOkeV); 

£n; : Energy of line photons (-8 keV) 

/?i : Internal radius of the shell, supposed to be com- 
parable to the stellar progenitor radius. In order 
to reduce the parameter space, the system will be 
solved for R t = 10 11 cm (typical WR radius) and 
R\ = 10 12 cm (typical blue supergiants radius). 

Li on : An efficiency of converting the illuminating con- 
tinuum into the sp ecific line ~ 1% is assumed 
(lLazzati et al.l20 02). For this reason the system will 
be solved for L; on ~ 10 49 erg s _1 ; 

/i e : Mean molecular weight. For a completely ionised 
material ju e = 1 , 1 .2 and 2 for pure hydrogen, solar 
composition and no hydrogen respectively. In the 
following we assume ji e « 2. 
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The unknown physical quantities of the system are: 



Funnel opening angle; 

Total Ni mass. W e expect MJ. < 1 M (see 
e.g. lArnettl Il996l) . In particular, if the unde- 
tected SN associated with GRB060904B is sim- 
ilar to other well studied SNe associated with 



GRB events, then 



0.4 M is ex- 



R 



pected (see e.g. Iwamoto et al. 1998; Mazza lFet alj 
l2003UMazzali et all 2006b: MazzalietalJ l2006al) . 
Results will be calculated for solar metallicity X = 
1.8xl0~ 3 ,X= 1.8 x 10~ 2 andX = 1.8 Xl0 _1 ; 
Outer radius of the shell. From the third relation of 
the equation [T4l we expect R - Ri> 10 12 cm; 
a : Mass density power-law index. The mass density 
of a SN shell can be approximated by a power- 
law in radius for radii greater than some critical 
value, while at radii less than critical the density is 
roughly constant. The critical radius evolves with 
time; its exact value is a function of the SN en- 
ergy, the ejected mass and the density power law 
index. However, at early times from the SN explo- 
sion (t < 1 hr) the density profile is nearly constant 
within R Q (see e.g. lBlondin & Ellisonl2001l) . At this 
dist ance the den sity undergoes a rapid cutoff ( see 
e.g. lArnettTl 996, their Fig. 13.2). For this reason we 
will consider only small values of the power-law in- 
dex: < a < 1. 

Our numerical results for the characteristic quantities of the 
problem are presented in Fig. [5H8] Figure [5] shows the varia- 
tion of the upper limit to the total mass of the SN shell derived 
from the second relation of the system ( TBI as a function of the 
funnel opening angle. As the temperature rises, the recombi- 
nation time increases, so that at a fixed funnel opening angle, 
more Ni is needed (and hence more mass in the SN shell) to 
fulfill the f rec = fj on requirement. Consequently the condition 
free > fion will be satisfied by a more extended range of SN 
shell masses: for example, for R\ = 10 n cm, R = 10 12 cm, 
T = 10 7 K we derive M < 1OOM with < 1.45 rad. For 
GRB-ass ociated HNe, the ejected mass reported in the literature 
(see e.g.llwamoto et alll998tlMazzali et alJl2003UMazzali et alj 
2006bt lMazzali et alj2006al) is however about one order of mag- 
nitude smaller. For this reason we expect M ~ 10 M . 

Figure [6] presents the variation of the total Ni mass required 
to explain the detected emission as a function of the funnel open- 
ing angle and the outer radius of the shell. Results are shown for 
different values of the internal radius and temperature. The strat- 
ification index has been fixed to 0. As the temperature rises, the 
required total Ni mass increases: this is due to the dependence 
of the recombination time on the temperature. On the other hand 
the increase of the funnel opening angle requires the detected 
emission to be explained by a smaller amount of Ni mass. This 
can be understood by considering that the density of the mat- 
ter in the SN shell is a monotonic increasing function of 0. 
Consequently f rec is a decreasing function of the funnel opening 
angle. The dependence of ML on the matter density manifests 
itself also through the parameter R a : the increase of the outer 
shell radius requires ML to increase, in order to reproduce the 
luminosity of the detected emission feature. Finally, increasing 
Ri, ML decreases. 

The dependence of the model on the value of the stratifica- 
tion index a is clear from Fig. [7] and [8] at a fixed outer radius, 
a — 1 requires smaller values of the funnel opening angle and 
higher values of the internal electron density. While the - R 



At 
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il-10keV 

(ergs -1 ) 



1 46 - 97 1.61^;; 8^ x 10 4!i 
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Wn 
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1 93 +0:n 

2 36+^ 
2 56 +0 - 12 
i a i +0.12 

2 3 f 5 M 

t 7^+8:1$ 
-Q.n 

9 99+0.15 



3 + o ? x 10 49 
2+° : i x 10 49 
2 + i\ x 10 49 

r°! x io 49 

4+^ x 10 48 
2^ x 10 45 



4^ 

-q. 

9+0 
z -0. 

2+o. 

r 

4 + 
2 + 



4 x 10 44 

\ x io 45 

I x 10 45 
J x 10 45 
1 x 10 45 
! x 10 45 
* x 10 44 



x 10 4 



Table 2. Rest frame isotropic 1 - lOkeV luminosity of 
GRB 060904B as a function of time from the GRB event. The 
XRT data have been splitted into different intervals of time, so 
that a minimum number of « 2000 photons is containd in each 
spectrum. A SPL model has been assumed. Af: interval of time; 
F: photon index (90% CI are listed); R: scale distance; «:particle 
density; ^ionization factor (see Sect. [7]). For the interval of time 
containing the emission feature (point 7) only the SPL contribu- 
tion has been reported. 



relation is nearly insensitive to the variation of this parameter, as 
a increases from to 1, the quantity n(X increases by a factor 
» 10 1 - 10 2 (Fig. [8]). Figure|7]shows the variations of the charac- 
teristic quantities of the problem assuming a total Ni mass of 0.4 
M : for T = 10 6 K, R Y = 10 11 cm, R = 10 12 cm the model re- 
quires 6 x 0.02 rad (a = 0). Increasing the electron temperature 
to 10 7 K (10 8 K) and keeping a = 0, the opening angle increases 
to 0.08 rad (0.48 rad). The corresponding internal electron den- 
sity at the emission site is n\ a 10 20 crrT 3 . The internal electron 
density would be about one order of magnitude lower if ju e ~ 16, 
value of the mean molecular weigth for a SN ejecta dominated 
by oxygen in a low ionization stage. 



5. The ionizing source 

The emission feature detected in GRB 060904B s pectra is char- 
acteris ed by a large EW — 2.0 + 0.3 keV. Following Lazzati et al. 
(2002), we expect an underlying 1-10 keV continuum luminos- 
ity Li-iokev > 100 Lii n e if this is to be identified with the ionizing 
source. From the Swift XRT light curve (see Fig. |2]i we derive 
^l-iokev = 4 x 10 48 ergs _1 at t ~ 160 s, time at which we de- 
tect the emission feature. The detected underlying continuum is 
therefore not physically linked to the line emission. 

The 1-10 keV rest frame luminosity as a function of time 
from the GRB onset is reported in Table [2] This table shows 
that the detected luminosity is actually higher than 100Lij ne ~ 
10 49 ergs _1 in the time interval 100-130 s. Assuming that the 
continuum illuminating the reprocessing material is similar to 
the observed radiation, it is therefore possible to explain the 
detected emission feature requiring the ionizing emission to be 
isotropic and the reprocessing material to have a large covering 
factor. In this scenario the time delay between the arrival of di- 
rect continuum photons and the observation of the claimed line 
emission A?d is explained in terms of light travel time effects. 

We assume that the ionizing photons are produced at t ~ 
100 s after the GRB onset or at a travelled distance of R ~ 3 x 
10 12 cm. At this distance they interact with the SN material. An 
estimate of the funnel opening angle 9 is therefore given by: 



(R + R p )e~ cAt D 



(15) 
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R I= 10 n (cm); T = 10 6 (K) R,= 10" (cm); T = 10 7 (K) R^IO 11 (cm); T = 10 8 (K) 




12 13 14 15 16 12 13 14 15 16 12 12.5 13 13.5 14 



Log (R„/cm) Log (R /cm) Log (R„/cm) 



Fig. 6. Total Ni mass contour plot according to the first relation of the system ([T4l > as a function of the funnel opening angle 9 and 
the outer radius of the shell R . Results are shown for the different electron temperatures indicated, R\ = 10 11 cm and 10 12 cm and 
a = 0. In each plot, from left to right, the black solid lines mark the M£ jj = KT 2 M , 10 _1 M Q , 1 M Q and 10 M G region. Red dashed 
lines: R = 10 12 cm and R - 10 13 cm. Lii ne = 1.10 x 10 47 erg s _1 , fNi = 8 keV have been used. 



R^IO 11 (cm) R I= 10 12 (cm) 




11.5 12 12.5 13 13.5 12.2 12.4 12.6 12.8 13 13.2 13.4 

Log(R Q /cm) Log(R /cm) 



Fig. 7. Funnel opening angle (6) versus outer shell radius (R ) according to Eq.[l4] first relation. A total Ni mass M^- = 0.4 M has 
been assumed (see the text for explanations). Results are shown for R\ = 10 11 cm (left panel) and R\ = 10 12 cm (right panel). In 
both panels the realtion is shown for three different temperatures: 10 6 K (bottom), 10 7 K and 10 8 K (top). Solid (dashed) black lines 
correspond to a — (a = 1). Dashed grey lines: R a = 10 12 , 10 13 cm. 



where R p is the progenitor radius. With A?d ~ 40 s, the previous 
relation translates into 6 ~ 22°, 17° and 5° for R p = 10 11 , 10 12 
and 10 13 cm. Moreover, assuming a typical expansion velocity 
Ksheii ~ 10 9 cms _1 , we have: 

Vshea(t + ?sn) ~ ct (16) 

where ?sn stands for the time delay between the GRB and the be- 
ginning of the SN expansion: for ?sn > the SN event happened 



before the GRB explosion. For t ~ 100 s the previous relation 
implies ?sn ~ 3000 s. 

6. Limits on the total GRB energy 

The emission feature detected in GRB 060904B shows an 
isotropic energy Lii ne X Af > 5 x 10 48 erg, where the in- 

equality accounts for the lower limit we have on the duration of 
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2. xlO 6 5.X 10 6 l.xlO 7 2.X 10 7 5. xlO'l.xlO 8 2. x 10 6 5. x 10 6 1. x 10 7 2. x 10 7 5.xl0 7 l.xl0 8 

T (K) T (K) 



Fig. 8. Internal electron density n t at the site of emission times the Ni abundance ratio X — M*/M as a function of the temperature. 
A total Ni mass ML = 0.4 M Q has been assumed (see the text for explanations). Left (right) panel: Ri = 10 11 cm (R t = 10 I2 cm). In 
both panels results are shown for different values of the outer radius of the SN shell R a . Solid (dashed) lines: a = (a = 1). 



the emission. Following Ghisellin i et al.1 (120021) the energy char- 
acterising the line emission can be used to set lower limits on 
the total energy reservoir of the GRB in the form of photons and 
kinetic energy: 

E =fZ> 200 " line (17) 

where rj x is the bolometric corretion while the symbol tj 7 stands 
for the efficiency in converting the kinetic energy into y-rays. A 
line efficiency of 1 % has been assumed. For GRB 060904B this 
translates into: 

E >_ 10"(*I)(M) erg (18) 

No amplification factor 4^/Qi; ne has been applied: from Se ct. [5] 
> 5°. Consequently, according to iGhisellini et all d2002l) the 
amplification factor due to electron scattering is ~ 1 . It is worth 
noting that the derivation of a numerical value for the r\ y pa- 
rameter, depends on many physical parameters of the burst itself 
and its environment: in the conventional internal+exte rnal shock 
GRB model r\ y ~ 1 % - 5% dZhang & Meszarosl2004 : however, 
allowing for an extre mely inhomogeneous velocity of the ejecta 
shell rj y ~ 40% (see iKobavashi et al.|[l997t iKobayashi & S aril 
1200 lb . In a more recent study IZhang et"aT] d2007al) derive for 
32 GRBs detected by Swift with early X-ray afterglow data 
rjy < 10% at the end of the shallow decay phase; r\ y > 90% 
at the deceleration time of the fireball. The bolometric correc- 
tion rj x is estimated using a Band spectrum dBand et alJI 1993b 
with typical values of the parameters: we obtain r\ x ~ 0.1. 

From equations [T7l and [T8l we derive E y ~ 10 52 (0. 1 /?7 r ) erg. 
An independent estimate of the isotropic equivalent radiated en- 
ergy Ey° can be obtained as follows. Using the correlation be- 
tween the SPL photon index F and the peak energy of GRB 
spectra measured by Swift-BAtT in the energy b and 15 - 150keV 
dZhang et alj|2007bi) the photon index of 1.70 dMarkwardt et alj 
120061) translates into E p ~ 85keV or E Pji ~ 15 0keV (intrin - 
sic peak energy). According to the Amati relation ( Amati 2006), 
the isotropic equivalent radiated energy in the 1-10000 keV cos- 
mological rest frame is therefore E y ° ~ 2 x 10 52 erg , in good 
agreement with the previous result. 

Moreover, assuming the time interval during which the ioniz- 
ing source has illuminated the reprocessing material to be com- 
parable with the duration of the detected emission (the recom- 
bination time would be too small to account for the duration 



of the emission), we have Em ~ L; on x At > 5 x 10 50 erg or 
Ej® T ~ 2 Li on x Af > 10 51 erg. The factor 2 in the previous re- 
lation has been added for geometrical reasons, since we expect 
a two-sided geometry with the line emitting material visible on 
only one side. For GRB 060904B this means that a fraction of 
the order of 1% of the radiated energy of the burst has illumi- 
nated the line emitting material. 



7. Discussion 

The XRT light curve of GRB 060904B is dominated by a giant 
flare in the first =s 200 s after the GRB onset (Fig. |2). At the 
end of this flare the spectrum becomes softer and an emission 
feature spectroscopically compatible with highly ionised Ni line 
emission arises. 

The detected emission can be explained in a reflection sce- 
nario (Sect. 0J. In these models high electron temperatures (> 
10 9 K) require total masses of potentially line emitting material 
as high as ~ 1OM for 9 < 60°. This would in turn imply a 
mass of the shell remnant > 1OOM even under the assumption 
of 100 times solar metallicity. For this reason, consistently with 
Sect. [3] we consider unlikely the T > 10 9 K scenario. If the fun- 
nel opening angle is to be identified with the jet opening angle 
Bj, the n ~ 60° at t ~ 10 2 s after the GRB onset is difficult to 
justify: Lamb et al. (2005), starting from theoretical considera- 
tions, concluded that GRBs have 9j ~ 1°; from the observational 
point of view, a typical 0, ~ 5° is derived from the afterglow 
jet break data dZhang et al.ll2004l) . If ~ 5° (~ 0.1 rad) the de- 
tected emission implies a total Ni mass of the order of 0. 1 M Q for 
Ri = 10 11 cm, R a = 10 12 cm and M* m > 1 M G for R a = 10 13 cm, 
with M N ;* increasing from T — 10 6 to 10 7 K. Such high val- 
ues of Ni mass suggest that the line emitting material is a SN 
remnant. (Note that up to this point no ad hoc assumption about 
the presence of the SN event has been made.) The same con- 
clusion is reached using Ri = 10 12 cm: in this case ~ 0.1 rad 
favours T ~ 10 6 K (see Fig. [6] and [7]) . We underline that in this 
scenario the giant flare that dominates the early XRT ligth curve 
can not be identified with the ionizing source: in Sect. [3] we de- 
rived > 5°. An undetected source of ionization is therefore 
required. Alternatively, non geometrical factors must be invoked 
to explain the delay between the arrival of direct continuum pho- 
tons and line photons. 

However, in most cases Swift data do not support the 
existence of standard jet breaks: using high quality multi- 
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wavelength data ICovino et alj d2006) concluded that no con- 
vincing case of achromatic break has been detected in Swift 
afterglows (the achromatic behaviour was not robustly es- 
tablished in the p ie-Swift era). T h e sam e conclusion has 
been reac hed by Panaitescu et al. d2006l) . Iwillingale et alj 
(2007j) and iBurrows & Racusinl ( 20071) in independent studies. 
Alternatively, the funnel opening angle might not be directly 
linked to the jet opening angle: in other words, the naive expec- 
tation 9 ~ 0j might be wrong. These reasons drive us to explore 
the prediction of the reflection model for larger 9. Moreover, in 
Sect. [5] we showed that the restarted activity of the GRB cen- 
tral engine might be the source of the observed line emission if 
the ionizing radiation is nearly isotropic and the material inter- 
cepts most of the ionizing flux. In that section we also showed 
that under the assumption that the reprocessing material is the 
SN shell expanding at v S h e n ~ 10 9 cms _1 , we obtain 9 ~ 22°, 
17° and 5° for progenitor radii of the order of 10 , 10 12 and 
10 13 cm, respectively. In this scenario the SN shell expansion 
began ~ 3000 s before the first GRB emission detected by BAT. 
Consequently, at t ~ 100 s we have R a ~ R p + 3 x 10 12 cm. 
With Mjjjj = 0.4 M Q , we find that the detected emission requires 
T ~ 10 7 K for 9 ~ 20° (0.3 rad) regardless of the internal ra- 
dius assumed. No prediction can be made about the stratification 
index. 

In this work we considered an efficiency of production of 
line photons by the reflection mechanism of the order of 1%. 
However, the line produ ction efficiency is a strong function of 
the ionization parameter dLazzati et alj|2002l) here defined as: 



4^D 2 F [ i_io] 



(19) 



where Fi_io is the 1 - lOkeV ionizing flux; Dl is the lumi- 
nosity distance (~ 4150 Mpc for z = 0.703); R is the distance 
of the reprocessing material and n is the particle density. The 
last column of Table [2] lists the quantity nR 2 for ^ = 10 4 ; the 
1 - lOkeV luminosity is derived from the XRT light curve. For 
a total shell mass of 10M G (order of magni tude of the stimated 
mass ejected by HNe explosion, see e.g. Iwamoto et al. 1998; 
IMazzali etaT] 120031; IMazzali et al.ll2006bllMazzali et al J 12006 a) 
and a scale distance R ~ 10 12 cm (~ v s h e ii X ?sn) we obtain 
nR 2 * 10 45 cm , value we actually derive from the observed 
XRT luminosity between 100 and 130 s for ^ = 10 4 . At this very 
high value of the ionization parameter, Ni line production effi- 
ciencies «J_%_ranbe obtained only for higher than solar metal- 
licities (Lazzati et al. 2002). In particular, the SN material asso- 
ciated with GRB 060904B must have a metallicity of the order 
of ten times the solar value. 

Finally, GRB jets might be structured, with angle-dependent 
energy per solid angle and Lorentz F factor as well: power-law 
jets or Gaussi an jets have been recent ly proposed in the liter- 
ature (see e.g. Zhang & Meszaros 2002). Unfortunately, the jet 
angular structure is still unknown: for this reason no assumption 
has been made about the angular distribution of the jet energy. 
As a consequence, our results have been calculated assuming 
that the detected radiation is representative of the whole emit- 
ted radiation, regardless of the direction of emission. However, 
it is important to underline that in the most realistic situation 
the emission properties of the ionizing continuum are likely to 
change with 9. 



8. Summary and Conclusions 

M07 recently reported the detection of a spectral feature 
in excess of the standard single power law behaviour in 



GRB 060904B spectra. The emission has been detected a 10 2 s 
after the onset of the GRB, lasted more than 47 s with an 
isotropic luminosity of * 10 47 erg s _1 (three orders of magnitude 
higher than the previous claimed detections). From the spectro- 
scopic point of view, the detected feature can be interpreted as 
He-like or H-like Ni emission. Given the poor statistics it is not 
possible to discuss the line profile. 

The peculiarity of the detection motivated us to explore 
in detail the theoretical implications of the presence of highly 
ionised Ni emitting =s 10 2 s after the burst. Given the typical 
tj » 1 (Sect. [3) only reflection models have been considered. 
Within these models, the external reflection scenario fails to pre- 
dict the observed luminosity when standard values of the wind 
velocity and mass loss rate of the progenitor star are used. On 
the other hand, this is Ni line emission, the detected feature 
can be interpreted in a funnel scenario with the line produced 
in reflection and typical Ni masses of the order of 0.1 M Q for 
9 = 9j ~ 5°, T ~ 10 6 K. Such high values of Ni masses strongly 
suggest the presence of a SN explosion associated with the GRB 
event. Alternatively, ;/ the reprocessing material is provided by 
the SN shell and ;/ the giant X-ray flare that dominates the early 
XRT light curve is the ionizing source, the expansion of the SN 
shell began « 3000 s before the GRB detection, M* ; ~ 0.4 M , 
9 ~ 20° and T ~ 10 7 K. These models imply mean electron den- 
sities n e ~ 10 21 cm 4 for ten times solar metallicity (X ~ 10" 2 ). 

In conclusion, starting from simple analytical models, we 
presented a consistent scenario in which the line emission can 
in principle arise. Given the importance of the claim, the subject 
certainly calls for further study. 

Acknowledgements. This work is supported by OAB-INAF by ASI grant 
1/011/07/0, by the Ministry of University and Research of Italy (PRIN 
2005025417) and by the University of Milano Bicocca (Italy). 



References 

Amati, L. 2006, MNRAS, 372, 233 

Arnett, D. 1996, in Supernova & Nucleosynthesis, ed. P. U. P. J. P. Ostriker 
Band, D., Matteson, J., Ford, L., et al. 1993, ApJ, 413, 281 
Blondin, J. M. & Ellison, D. C. 2001, ApJ, 560, 244 
Bottcher, M. 2004, Advances in Space Research, 34, 2696 
Burrows, D. N. & Racusin, J. 2007, ArXiv Astrophysics e-prints 
Butler, N. R. 2007, ApJ, 656, 1001 

Campana, S., Mangano, V., Blustin, A. J., et al. 2006, Nature, 442, 1008 
Covino, S., Malesani, D., Tagliaferri, G., et al. 2006, ArXiv Astrophysics e-prints 
Dickey, J. M. & Lockman, F. J. 1990, ARA&A, 28, 215 

Fugazza, D., D' Avanzo, P., Malesani, D., et al. 2006, GRB Coordinates Network, 
5513, 1 

Galama, T. J., Vreeswijk, P. M., van Paradijs, J., et al. 1998, Nature, 395, 670 

Ghisellini, G., Lazzati, D., Rossi, E., & Rees, M. J. 2002, A&A, 389, L33 

Iwamoto, K., Mazzali, P. A., Nomoto, K., et al. 1998, Nature, 395, 672 

Kobayashi, S., Piran, T., & Sari, R. 1997, ApJ, 490, 92 

Kobayashi, S. & Sari, R. 2001, ApJ, 551, 934 

Lamb, D. Q., Donaghy, T. Q., & Graziani, C. 2005, ApJ, 620, 355 

Lazzati, D. 2000, Phd thesis, 999 

Lazzati, D., Ramirez-Ruiz, E., & Rees, M. J. 2002, ApJ, 572, L57 
Malesani, D., Tagliaferri, G., Chincarini, G., et al. 2004, ApJ, 609, L5 
Markwardt, C., Barbier, L., Barthelmy, S. D., et al. 2006, GRB Coordinates 

Network, 5520, 1 
Mazzali, P. A., Deng, J., Nomoto, K., et al. 2006a, Nature, 442, 1018 
Mazzali, P. A., Deng, J., Pian, E., et al. 2006b, ApJ, 645, 1323 
Mazzali, P. A., Deng, J., Tominaga, N, et al. 2003, ApJ, 599, L95 
Moretti, A., Margutti, R., Pasotti, E, et al. 2007, ArXiv e-prints, 711 
O'Brien, P. T„ Willingale, R., Osborne, J., et al. 2006, ApJ, 647, 1213 
Panaitescu, A., Meszaros, P., Burrows, D., et al. 2006, MNRAS, 369, 2059 
Patat, F, Cappellaro, E., Danziger, J., et al. 2001, ApJ, 555, 900 
Piran, T. 2005, Reviews of Modern Physics, 76, 1143 

Protassov, R., van Dyk, D. A., Connors, A., Kashyap, V. L., & Siemiginowska, 

A. 2002, ApJ, 571,545 
Sako, M., Harrison, F. A., & Rutledge, R. E. 2005, ApJ, 623, 973 
Stanek, K. Z., Matheson, T., Garnavich, P. M., et al. 2003, ApJ, 591, L17 



10 



R. Margutti et al.: Anomalous X-Ray emission in GRB 060904B: a Nickel line? 



Sunyaev, R. A. & Titarchuk, L. G. 1980, A&A, 86, 121 

Verner, D. A. & Ferland, G. J. 1996, ApJS, 103, 467 

Verner, D. A. & Yakovlev, D. G. 1995, A&AS, 109, 125 

Vietri, M., Ghisellini, G., Lazzati, D., Fiore, F., & Stella, L. 2001, ApJ, 550, L43 

Willingale, R., O'Brien, P. X, Osborne, J. P., et al. 2007, ApJ, 662, 1093 

Woosley, S. E. & Weaver, T. A. 1995, ApJS, 101, 181 

Zhang, B., Dai, X., Lloyd-Ronning, N. M., & Meszaros, P. 2004, ApJ, 601, LI 19 
Zhang, B., Liang, E., Page, K. L., et al. 2007a, ApJ, 655, 989 
Zhang, B. & Meszaros, P. 2002, ApJ, 571, 876 

Zhang, B. & Meszaros, P. 2004, International Journal of Modern Physics A, 19, 
2385 

Zhang, B., Zhang, B.-B., Liang, E.-W., et al. 2007b, ApJ, 655, L25 



